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ABSTRACT 

We discuss the effect of dispersion delay due to intervening ionized media in the radio 
emission of gamma-ray bursts (GRBs). For high redshift GRBs (z <; 3), the ionized 
intergalactic medium (IGM) should be the dominant source of dispersion without 
substantial local or foreground contamination, offering a unique probe of the cosmic 
reionization epoch through measures of the free electron column density out to different 
redshifts. The delay times from z ~ 10 can be ~ 1 hour at 100 MHz and ~ 10 hours 
at 30 MHz. On the other hand, dispersion by local ionized material may be important 
for GRBs at lower redshifts if they occur inside or behind dense molecular clouds, 
providing clues to the GRB environment; free- free absorption may also be significant 
in this case. Detecting dispersion delay in the known radio afterglow emission should 
be extremely challenging due to the low fluxes at the relevant frequencies, but may 
be marginally possible for rare, bright afterglows by the Square Kilometer Array. If 
GRBs also emit prompt, coherent radio emission as predicted by Sagiv & Waxman 
(2002), the observational prospects can be much better; detection of a sufficient sample 
of high z GRBs by the Low Frequency Array may allow a discrimination of different 
reionization histories. Interesting constraints on the low redshift, warm-hot IGM may 
also be obtainable through dispersion. 

Key words: cosmology: theory - intergalactic medium - gamma-rays: bursts - ra- 
diation mechanisms: non-thermal - radio continuum: general 



1 INTRODUCTION 

Gamma-ray bursts (GRBs) are being increasingly recog- 
nized as potentially powerful probes of the very high redshift 
universe (Loeb 2003, Djorgovski et al. 2003). The mount- 
ing evidence that GRBs (at least with long durations) are 
closely associated with the formation and collapse of massive 
stars (e.g. van Paradijs, Kouveliotou & Wijers 2000, Bloom, 
Kulkarni & Djorgovski 2002, Hjorth et al. 2003, Matheson et 
al. 2003) indicates that they can occur in significant numbers 
at early epochs, perhaps up to the formation era of the very 
first stars (Lamb & Reichart 2000, Bromm & Loeb 2002). 
Their bright, broadband emission, a large part of which can 
be robustly modeled as nonthermal emission from a relativis- 
tic blast wave, should be readily detectable out to the high- 
est redshifts they are expected to exist (Lamb & Reichart 
2000, Ciardi & Loeb 2000). Given the likelihood that bright 
quasars are much rarer at high z, GRBs may play unique 



* inouemu@MPA-Garching.MPG.DE 
© 2003 RAS 



roles as beacons illuminating the so-called "dark ages" of 
the universe (Barkana & Loeb 2001, Miralda-Escude 2003). 

An issue of considerable current interest is the epoch of 
cosmic reionization. The question of when and how reioniza- 
tion in the universe occurred and what caused it is critical, 
having profound implications for the formation of the first 
stars and galaxies and their subsequent evolution (Madau 
& Kuhlen 2003, Haiman 2003). The highest redshift quasars 
discovered by the Sloan Digital Sky Survey (SDSS) have re- 
vealed complete spectral absorption troughs (Gunn & Pe- 
terson 1965), showing that the neutral fraction of the in- 
tergalactic medium (IGM) is at least ~ 10~ 3 at z ~ 6 
(Djorgovski et al. 2001, Becker et al. 2001, White et al. 
2003). On the other hand, measurements of temperature- 
polarization correlations in cosmic microwave background 
anisotropies by the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite indicate that the Thomson optical depth 
to the last scattering surface is r e ~ 0.16, which in turn is 
suggestive of a high reionization redshift z ~ 20 (Kogut et al. 
2003, Spergel et al. 2003). Various scenarios have been put 
forth to account for these observations which predict very 
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different histories for the evolution of the ionization frac- 
tion (e.g. Ciardi, Ferrara & White 2003, Cen 2003, Wyithe 
& Loeb 2003; see Miralda-Escude 2003, Haiman 2003 for 
reviews). GRBs may be the most promising tools to differ- 
entiate between these possibilities, and some methods using 
near-infrared spectroscopy (Miralda-Escude 1998, Barkana 
& Loeb 2003) or photometry (Inoue, Yamazaki & Nakamura 
2003) have been proposed. 

On an entirely different note, the properties of the lo- 
cal environment in which GRBs explode is a crucial issue, 
bearing on the still mysterious nature of the GRB progen- 
itors. At present, the density of the medium surrounding 
GRBs is quite uncertain, even for the well-observed bursts 
at low redshifts. Through self-consistent broadband mod- 
eling of a sample of bright GRB afterglows, Panaitescu & 
Kumar (2001, 2002) derive ambient densities in the range 
10~ 3 - 30cm" 3 on the scales of external shock decel- 
eration (radii r ~ 10 16 — 10 18 cm). In light of the numerous 
indications for massive star progenitors, one may expect a 
denser, radially decreasing medium characteristic of stellar 
winds (Chevalier & Li 1999), but this is generally not sup- 
ported by observations. On the other hand, there is tentative 
(but far from conclusive) evidence from soft X-ray absorp- 
tion and optical/X-ray flux ratios that the GRB environ- 
ment may have higher densities characteristic of Galactic 
molecular clouds (MCs), n ~ 10 2 — 10 3 cm -3 or even larger 
(Galama & Wijers 2000, Reichart & Price 2002, Piro 2002 
and references therein). The actual environment may have 
a complex density structure such as that of superbubbles 
(Scalo & Wheeler 2001). Clearly, one would like to have as 
many ways as possible to probe the nature of the circum- 
burst medium. 

We make the case here that dispersion delay in the ra- 
dio emission of GRBs may be capable of providing valuable 
insight into both of these issues, t Dispersion is sensitive to 
the integrated column density of ionized material along the 
line of sight, and is therefore unique and complementary to 
methods utilizing absorption by intervening neutral gas. We 
show that dispersion delay may be effective in discriminating 
between different reionization histories, particularly if GRBs 
emit strong, coherent low frequency emission as predicted 
by some recent models (Sagiv & Waxman 2002). Interest- 
ing constraints may also be achievable for the ill-understood 
warm-hot IGM at low redshift (Cen & Ostriker 1999, Dave 
ct al. 1999, 2001). Dispersion of GRB radio emission as a 
probe of the IGM was first discussed by Ginzburg (1973) 
and later by Palmer (1993). We adopt the following cosmo- 
logical parameters that are consistent with WMAP results: 
n m = 0.27, n A = 0.73, h = 0.71 and n b = 0.044 (Spergel et 
al. 2003). 

This paper is organized as follows. In §2, we evaluate 
the degree of dispersion due to the ionized IGM and the dif- 
ferences expected for various proposed reionization scenar- 
ios. Dispersion due to different types of local ionized media 
and its evolution with redshift is discussed in §3. §4 ad- 
dresses free-free absorption in these environments. In §5, we 
discuss the observability of dispersion delay in the known 
afterglow emission from GRBs as well as the predicted co- 
herent emission. We conclude in §6 along with remarks on 

t A similar study has recently been carried out by loka (2003). 



future prospects. The formalism for calculating the radio 
afterglow emission is detailed in the Appendix. 



2 DISPERSION IN THE IONIZED 
INTERGALACTIC MEDIUM 

An electromagnetic signal of frequency v will propagate 
through an ionized medium with a velocity v = c(l — 
Vp /v 2 ) 1 / 2 , less than the speed of light in vacuum (e.g. Ry- 
bicki & Lightman 1979). Here v v = (n e e 2 /7rm e ) 1/ ' 2 is the 
plasma frequency and n e is the free electron density of the 
medium. Therefore the arrival time of a low frequency sig- 
nal will be delayed with respect to that at v S> v p by a 
frequency-dependent amount which is proportional to the 
integrated column density of free electrons along the propa- 
gation path, 

where dl is an element of distance along the path. The 
dispersion measure (DM) J n e dl is commonly expressed in 
units of pc cm~ 3 . Observations of dispersion in the radio 
emission of pulsars in our Galaxy show that the DM of the 
Galactic disk is ~ lOOpc cm~ 3 for typical lines of sight, with 
maximum values of ~ lOOOpc cm -3 for lines of sight near 
the Galactic center (e.g. Cordes & Lazio 2002). 

The intergalactic medium (IGM) is known to be sig- 
nificantly ionized out to at least z Ss 6 (Gunn & Peterson 
1965, Djorgovski et al. 2001). The radio emission from high 
redshift GRBs may then be affected by dispersion in the 
ionized IGM (Ginzburg 1973, Palmer 1993). Note that the 
discussion of this section may also apply to any other type 
of transient radio source. For propagation through cosmo- 
logical distances, the appropiate comoving distance element 
is dl = c\dt/dz\dz, where \dt/dz\ = (1 + z)~ H(z)^ 1 and 
H{z) = (fi m (l + zf + Ov) 1/2 for a fiat cosmology, t One 
must also account for the redshifting frequency and time di- 
lation during propagation. Thus, for a source at redshift z, 
the mean dispersion delay at observer frequency v due to 
the IGM is 

A+ / \ e 2 1 f z , dt n e (z') . . 

AtlGM( ^' 2) = J dZC d?—- (2) 

The mean free electron density is n e (z) = x e {z)n e fl{l + z)' i , 
where x e (z) is the ionization fraction. We limit the discus- 
sion here to the average effects of the ionized IGM assum- 
ing homogeneity; the implications of a more realistic, inho- 
mogeneous IGM structure is touched upon in §6. Assum- 
ing that the IGM comprises Y = 0.24 of helium by mass 
fraction and the rest hydrogen, and that x e (z) of the hy- 
drogen and helium are fully and singly ionized, respectively 
(e.g. Yoshida et al. 2003), n e , = 0.92 x lO^n^cm -3 , or 
n e , = 2.1xl0~ 7 cm~ 3 using the WMAP value of Q b h 2 . (He- 
lium may become fully ionized at z 3, but this is neglected 
as it will only result in a small difference.) 

The DM of the IGM out to different redshifts, defined 

t Strictly speaking, the propagation path will deviate from a null 
geodesic, but this effect is negligibly small unless v <C c. 
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appropriately as the integral in Eq.2, is shown in Fig. la as- 
suming for illustration that x e (z) = 1, i.e. the universe is al- 
ways ionized. We see that the IGM DM already exceeds typ- 
ical Galactic values at redshifts as low as z ~ 0.2, surpasses 
~ lOOOpc cm -3 at z ^ 1, and can reach ~ 10 4 pc cm -3 at 
2 J; 16. Fig. lb shows the dispersion delay times at different 
observer frequencies. At v =300, 100 and 30 MHz, the cor- 
respsonding delay times for sources at z ~ 10 are roughly 5 
minutes, 1 hour and 10 hours, respectively. Thus, dispersion 
delay may cause a potentially measurable imprint of propa- 
gation through the ionized IGM in the light curves of GRB 
radio emission, as we discuss in §5. 

We also see that the dependence of DM and At with 
respect to z becomes much milder for z J> 5. This results 
from the v~ 2 dependence of dispersion and the associated 
(1 + z)~ 2 factor in the integral of Eq.2, which counteracts 
the (1 + z) 3 increase in the density and (1 + z) increase 
in time dilation: at a given observer frequency, the higher 
the source redshift, the higher the rest-frame frequency, and 
hence less dispersion. This tendency is not so favorable for 
the purpose of probing the details of the reionization epoch, 
for which the most interesting range of redshifts are z ~ 
5 — 20. Nevertheless, different possible reionization histories 
can lead to observationally significant differences in At and 
its z-dependence, which may give us a crucial handle on 
discriminating between various scenarios for reionization. 

To demonstrate this, we consider two different types 
of reionization scenarios which have been discussed in light 
of the WMAP results. First, we take a more conventional 
picture of reionization by early massive stars, referring 
to the numerical simulation results of Ciardi, Ferrara & 
White (2003) for concreteness (see also Miralda-Escude 
2003, Haiman 2003 and references therein). Making rela- 
tively conservative assumptions about the properties of the 
ionizing stars and their initial mass function (IMF), these 
models predict that the universe is gradually ionized un- 
til an epoch of complete reionization z r is reached, the 
value of which depends on the form of the IMF, produc- 
tion efficiency of ionizing photons, and their escape frac- 
tion from the host galaxy. The numerically calculated evo- 
lution of the ionization fraction is given in Fig. 3 of Ciardi, 
Ferrara & White (2003). We schematically describe this as 
x e (z) = io-°- 3 C*-*r) f or z > Zr an d Xe ( z ) = 1 for z < z r , 
and take the cases of z r =6, 9 and 13.5. The latter two cases 
correspond to their models S5 and L20, which are consistent 
with the WMAP Thomson optical depth but not with the 
z <; 6 quasar Gunn-Peterson troughs. On the other hand, 
the additional case of z r =6 considered here, which may be 
realized with a lower photon production efficiency, may be 
more concondant with the z ^ 6 quasar observations, but 
less with the WMAP results. 

Second, we consider the scenario proposed by Cen 
(2003) (see also Wyithe & Loeb 2003, Sokasian et al. 2003), 
in which the universe is reionized twice: first by very mas- 
sive Population III stars, after which a period of partial 
recombination follows, and then for the second time by 
Population II stars with a different IMF and photon pro- 
duction efficiency. Fig. 9 of Cen (2003) shows the ioniza- 
tion fraction in his model, which we schematically represent 
as x e (z) = io-°- 3 ( z - 16 ) for z > z\ = 16, x e (z) = 1 for 
13 < z < 16, x e {z) = 10" ' 3 for 6 < z < 13 and x e {z) = 1 
for z < Z2 — 6. Note that in the intermediate recombination 
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Figure 1. a) Dispersion measures for dilferent source redshifts 
due to the IGM (solid curve), molecular cloud environments of 
n MC = 10 2 cm — 3 (lower dashed curve) and n^c = 10 3 cm -3 
(upper dashed curve), and the host galaxy ISM with tijsm = 
(1 + z) 3 cm -3 (dotted curve). The curves ending horizontally at 
high z correspond to different reionization histories as discussed in 
the text and Fig. 2. The typical Galactic dispersion measure (dot- 
dashed line) is also shown, b) Dispersion delay times in hours due 
to the IGM (solid curves) and MCs of n M c = 10 3 cm -3 (dotted 
curves), for observer frequencies v =10, 30, 100 and 300 MHz, 
from top to bottom. 



period, the ionization fraction is still considerable, but the 
amount of neutral hydrogen is enough to make the universe 
highly opaque to Lya photons (Cen 2003). 

The dispersion delay times at i/=30 MHz for these dif- 
ferent reionization histories are displayed in Fig. 2 (note the 
linear vertical scale). For the gradual reionization scenarios, 
At keeps increasing with z until the reionization redshift z r 
is reached, and remains constant thereafter at z z r . For 
the two-epoch reionization scenario, the trend of At with 
z becomes shallower above the second reionization redshift 
but does not completely level off as the recombination in 
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Figure 2. The dispersion delay times in hours due to the IGM 
at v =30 MHz for scenarios with gradual reionization (dotted 
curves) and two epoch reionization (short dashed curve). The 
case for a perpetually ionized IGM (solid curve) is also shown. 



the intermediate epoch is only partial. The behavior steep- 
ens again from the second recombination redshift until the 
first reionization redshift, above which At becomes constant. 

For ^=30 MHz, At at z ~ 10-15 is ~ 10 hours, and the 
the different reionization histories lead to differences in At of 
several hours. At z/=100 MHz, these are ~ 1 hour and tens 
of minutes, respectively. If the radio emission of GRBs at 
these frequencies are bright enough, these delays ought to be 
measurable (see §5.2). Measuring At for a number of GRBs 
occurring at different redshifts, and determining their red- 
shifts independently, e.g. through infrared Lya breaks (Cia- 
rdi & Loeb 2000, Gou et al. 2003), X-ray lines (Meszaros & 
Rees 2003) or 21 cm absorption features (Furlanetto & Loeb 

2003) § in their afterglow spectra, we may be able to effec- 
tively discriminate between different reionization histories. 
GRB radio emission may then serve as extremely valuable 
probes of the cosmic reionization epoch which is otherwise 
difficult to investigate directly. 

We mention that below z ~ 3, a progressively larger 
fraction of the IGM mass should be incorporated into the 
warm-hot phase at temperatures T ~ 10 5 — 10 7 K, being 
shock heated by infall onto incompletely virialized large scale 
structure such as filaments (Cen & Ostriker 1999, Dave et al. 
1999, 2001). Whereas the average overdensity of the cooler, 
photoionized IGM giving rise to the Lya forest should be 
8 ~ 0, this warm-hot IGM (WHIGM) is expected to possess 
overdensities S ~ 10 — 30, systematically above the mean 
baryon density of the universe (Dave et al. 2001). Although 
dispersion is insensitive to the medium temperature, this 
difference in the average density of the WHIGM may cause 



3 However, we find that radio afterglows from high-z bursts are 
typically too faint for detecting 21 cm absorption features due to 
the host galaxy; §5.1. 



an increase in the dispersion above the simple prediction of 
Eq.2, especially as the contribution to the DM from z 3 is 
very large (Fig.l). However, even though the mass fraction 
of the WHIGM can be as high as 30-40 % at z ~ (5-20 % 
at z ~ 2) (Dave et al. 2001), its volume fraction is expected 
to be much smaller, less than 10 % at z ~ (5 % at z ~ 2) 
(Cen & Ostriker 1999). Thus, the fraction of the line of sight 
to a GRB that is occupied by the higher density WHIGM 
may not be so large. Nevertheless, this contribution from the 
WHIGM may be important for some lines of sight, possibly 
complicating studies of the reionization epoch with high-z 
bursts. On the other hand, given the difficulty of directly ob- 
serving this important but poorly understood component of 
the low redshift IGM (compare e.g. Tripp, Savage & Jenk- 
ins 2000, Mathur, Weinberg & Chen 2003, Zappacosta et 
al. 2002), interesting constraints may be obtainable through 
dispersion in GRBs at z Ss 3 (provided that local dispersion 
is unimportant; §3), particularly if the degree of inhomo- 
geneity can be quantified (see §6). Further discussion of the 
effect of the WHIGM on dispersion delay is beyond the scope 
of this paper, but may warrant a more detailed study. 



3 DISPERSION IN IONIZED LOCAL 
ENVIRONMENTS 

Besides the global IGM, there may be ionized regions lo- 
cally enshrouding the GRB that can give rise to additional 
dispersion. Such local contributions may be described by 



At\ oc (v, z) = 



1 N e (z) 



2irm e c v 2 1 + z 



(3) 



where N e (z), the column density of free electrons local to 
the source, can in general evolve with z. The DM here is 
N e (z)/(l + z). 

Local ionized regions can come about in two ways: i) 
photoionization of the GRB vicinity by UV photons from 
the GRB itself (e.g. Perna & Loeb 1998, Fruchter, Krolik & 
Rhoads 2001, Draine & Hao 2002, Perna & Lazzati 2002), or 
ii) photoionization (and/or mechanical heating) of the GRB 
surroundings by a population of massive stars preexisting 
in the host galaxy (e.g. Ciardi et al. 2001, Barkana & Loeb 
2003). 

In the second case, the HII regions may extend well out 
into the ambient IGM for galaxies above the reionization 
redshift. However, we can immediately see that the disper- 
sion from the local HII regions in the IGM should be gener- 
ally miniscule. Depending on the IMF and age of the massive 
star population, the radii of these regions should be in the 
range run ~ 100 — 300 kpc (Barkana & Loeb 2003). Since 
the IGM density is n IG M ^ 2.1 x 10~ 7 cm~ 3 (l + z) B , the 
DM is n IG Mr H ii/(l + z) ~ (2 - 6) x 10~ 2 pc cm~ 3 (l + zf , 
below the Galactic DM even at z ~ 30. As one approaches 
the global reionization epoch, the HII regions become much 
larger and begin to mutually overlap, but then the associ- 
ated dispersion should be attributed to the global IGM. 

Local dispersion can be important for much denser loca- 
tions, as discussed below. In such regions, preexisting mas- 
sive stars play a lesser role compared to GRB UV emission 
in photoionizing the environment, since the hydrogen recom- 
bination time tree ~ 10 5 yr n _1 can be much shorter than 
the lifetime of the massive stars, 10 7 — 10 s yr. 
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3.1 Molecular Clouds 

The current controversy regarding the environments of the 
observed GRBs and the circumstantial evidence for molec- 
ular clouds was described in §1. As an exemplary case, we 
here consider the dispersion that may be caused by environ- 
ments similar to Galactic MCs, with typical size tmc ~ 10 
pc and densities umc ~ 10 2 — 10 3 cm~ 3 (e.g. Reichart & 
Price 2002; note that MC cores may have smaller sizes and 
higher densities) . Under the assumption that the GRB emits 
a flash of 10 50 £W,50 erg in the UV (Draine & Hao 2002; al- 
though the actual UV fluence may vary considerably from 
burst to burst; see remark in §4), these UV photons can 
ionize hydrogen with density n out to a radius of roughly 

Uon ~ 10 20 cm Bj/v iB0 n" 1/3 . (4) 

For n M c ~ 10 2 - 10 3 cm~ 3 , r ion ~ 3 - 7 pc. If the MC 
properties do not evolve with z, the corresponding DMs will 
be ~ (720 - 3300) (1 + z) _1 pc cm~ 3 . Compared with the 
IGM DM and delay times in Fig.l, we see that dispersion 
due to MCs may dominate at low redshifts z <, 1 — 2, but 
becomes rapidly unimportant at higher z, simply due to the 
(1 + z)- 1 factor. Thus, by measuring this effect for low- z 
GRBs one may gain valuable insight into the controversial 
issue of the GRB environment (although free- free absorption 
may foil this prospect; see §4). At the same time, such local 
contributions will not contaminate dispersion by the IGM 
in high-z GRBs. 

3.2 Interstellar Medium of Host Galaxy 

At high redshifts, the deficiency of metals and the associ- 
ated reduction in cooling may thwart the formation of MCs 
similar to the Galactic ones. On the other hand, the average 
density of the interstellar medium (ISM) of host galaxies 
should increase with z, possibly leading to another source of 
dispersion. Lacking observational constraints, there is con- 
siderable uncertainty in how the ISM density evolves at high 
z. The properties may be different for the metal-free Popu- 
lation III hosts and the more metal-enriched Population II 
hosts, the transition between the two possibly taking place 
around z ~ 10 - 15 (Cen 2003, Wyithe & Loeb 2003, Sal- 
vaterra & Ferrara 2003). 

Based on the disk galaxy models of Fall & Efstathiou 
(1980) and Mo, Mao & White (1998), which are consis- 
tent with the observed properties of disks and damped Lya 
absorbers at z Si 3, Ciardi & Loeb (2000) have discussed 
how the density of high z disk galaxies may evolve. In 
their model, the average ISM density and the disk scale 
height evolve respectively as uism oc M 2 ^ 3 (l + z) A and 
h tx M _1 ' 3 (l + z)~ 2 , where M is the mass of the galaxy 
(which is proportional to that of the dark matter halo). 
Thus, for a fixed galaxy mass, the disk becomes very dense 
and thin at high z. However, the standard picture of hier- 
archical structure formation implies that the typical mass 
of galaxies should also evolve strongly. For example, taking 
the minimum star-forming halo mass to be M m j n = 4.4 x 
1O 9 M (1 + z)- 1 - 6 ^- 1 (Ciardi & Loeb 2000), the halo mass 
averaged over the Press-Schechter mass function evolves 
from M - 1.4 x 10 n M Q at z = to M ~ 4.7 x 1O 8 M 
at z = 10, decreasing by more than two orders of magni- 
tude. The actual mass range relevant for GRB hosts should 



be that weighted by the amount of star formation, the de- 
tails of which is uncertain and model-dependent (e.g. Ciardi 
& Loeb 2000, Barkana & Loeb 2003). In any case, if one 
is concerned with the typical GRB interstellar environment, 
the density evolution in this model should be much milder 
than simply taking (1 + z) 4 (c.f. Gou et al. 2003). 

On the other hand, if the disk forms by collapsing by a 
constant factor in density with respect to that of the halo, 
the ISM density should scale as (1 + z) 3 , as argued by Cen 
(2003). The dependence on mass through the typical halo 
collapse redshift may be weak. If we follow Cen (2003) in 
normalizing to the Galactic ISM density, the density evolu- 
tion in this case is uism — (1 + z) cm . This may possi- 
bly be consistent with numerical simulation results of pri- 
mordial galaxy formation which sometimes reveal disk-like 
structures with densities n ~ 10 3 — 10 4 cm -3 at z ~ 20 
(Abel, Bryan & Norman 2002, Bromm, Coppi & Larson 
2002, Yoshida et al. 2003). 

In either situation, the DM due to the host galaxy 
ISM should be subdominant. For the (1 + z) 3 case, the ion- 
ization radius due to the GRB photons would be r ion ~ 
33pc(l + z) _1 , so that the DM is ri 0n nisM(z)/(l + z) ~ 
33pc cm~ 3 (l + z), well below the IGM DM (Fig.l). Even 
if one makes the extreme (and unrealistic) assumption of 
nisM = (l + z) 4 cm~ 3 , the DM is 33pc cm~ 3 (l + z) 5/3 , which 
can reach 10 4 pc cm~ 3 at z ~ 30 but is always overwhelmed 
by the IGM dispersion. 

In addition, radiative and/or mechanical feedback, ei- 
ther by the GRB progenitor star or preexisting massive 
stars, may be strong enough to keep the ambient densities as 
low as n ~ 10~ 2 - lcm~ 3 (Scalo & Wheeler 2001, Whalen, 
Abel & Norman 2003), regardless of the average ISM density 
of the host galaxy, further diminishing the DM. 

Summarizing up to this point, molecular cloud envi- 
ronments can be the main contributor to dispersion at low 
redshifts z ^ 1 — 2, while the mean IGM should dominate 
at higher redshifts. Other possible sources of dispersion are 
all expected to be minor. 



4 FREE-FREE ABSORPTION 

Before turning to the observability of radio dispersion in 
GRBs, we discuss the effect of free-free absorption, which 
can strongly affect low frequency radio emission in dense, 
ionized environments, as is the case for thermal emission 
from Galactic HII regions. For an ionized medium consisting 
of pure hydrogen, density ncm -3 , temperature IQ^T^K and 
depth r pc pc, the optical depth to free-free absorption at 
observer frequency v = 10 8 i/8 Hz is 

ty, = 5.5 x 10~ V 2 (l + z)~ 2 n 2 T- 3/2 r pc gff , (5) 

where the Gaunt factor appropriate for the radio regime is 
g ff = 0.281n(T 4 ///f(l+z) 2 )-7.2 (Rybicki & Lightman 1979, 
Shu 1991). 

Free-free absorption by the IGM is completely negligi- 
ble, as t„ ~ lO -8 ^ 2 even when assuming a constant tem- 
perature of Ti = 1; in reality the temperature should rise 
at higher z (e.g. Carilli, Gnedin & Owen 2002). However, 
absorption can be serious for the dense local environments. 
For the MCs discussed above with umc = 10 2 — 10 3 cm -3 
and rion ~ 3 — 7 pc, r v ~ 0.7 — 30i^~ 2 at z = 1. For the host 
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galaxy ISM with uism = (l+z) 3 cm 3 , r ion ~ 33pc(l+.z) 1 , 
so that t v = 1.4 x 10" 3 i/ g ~ 2 (l + z) 3 T~ 3/2 g ff , which is 
n, ~ lAvg 2 at z = 10 and t„ ~ 30^ 2 at z — 30. Un- 
der such circumstances, radio emission below 0.1-1 GHz can 
be strongly attenuated, hampering the chances of observing 
dispersion effects at lower frequencies. On the other hand, 
observing the spectral break due to free-free absorption may 
be used to advantage as an alternative probe of dense, ion- 
ized environments. 

One loophole is that the efficiency of flash photoioniza- 
tion of the GRB surroundings may vary significantly from 
burst to burst, as the fluence of UV ionizing photons emit- 
ted from the reverse shock can depend strongly on the ini- 
tial Lorentz factor, GRB duration and other parameters 
(Kobayashi 2000, Waxman & Draine 2000, Draine & Hao 
2002, Perna & Lazzati 2002). In the case of inefficient pho- 
toionization by the GRB, neither free-free absorption nor 
dispersion by the local environment will affect the radio 
emission, unless the prior activity of nearby massive stars 
have been vigorous. Dispersion by the IGM should be the 
only effect of observational consequence. 



5 DISPERSION DELAY IN GRB RADIO 
EMISSION 

We now discuss the observability of dispersion delay effects 
in the radio emission from GRBs. 

5.1 Radio Afterglows 

We first consider the known radio afterglow emission, which 
can be well-described as incoherent synchrotron emission 
from electrons accelerated by a relativistic blastwave decel- 
erating in the ambient medium (Paczyhski & Rhoads 1993, 
Katz 1994, Meszaros & Rees 1997, Sari, Piran & Narayan 
1998; see Piran 1999, van Paradijs et al. 2000, Meszaros 2002 
for reviews) . The calculation of the spectra and light curves 
follows standard discussions and is detailed in the Appendix. 
Particular attention has been made to account correctly for 
synchrotron self-absorption at high ambient densities. Al- 
though radio emission is expected from both the long-lived 
forward shock and the shorter-lived reverse shock (Sari & 
Piran 1999, Kobayashi & Sari 2000), we neglect the latter 
component as we find it to be important only above a few 
GHz. At the lower frequencies of observational interest for 
dispersion effects, the forward shock emission is always dom- 
inant. 

The afterglow flux at radio frequencies depends strongly 
on the density of the ambient medium, especially in the self- 
absorbed part of the spectrum. Much of the discussion of §3 
regarding the density of the GRB environment and its high 
z evolution is also relevant here. However, we note that the 
ambient density on the scales of external shock deceleration 
(radii 10 16 - 10 18 cm) may or may not be different from that 
contributing to the dispersion (and free- free absorption) . For 
example, the external shock may happen to decelerate in a 
location of low density while the afterglow emission passes 
through a much denser region in the foreground toward the 
observer. This may occur in a complicated geometry like 
that of superbubbles (Scalo & Wheeler 2001), or even for 
Population III conditions if feedback from the progenitor 



V=100 MHz, t=0.01 day- 




z 



Figure 3. The dependence on z of the afterglow flux at v = 
100 MHz and t = 0.01 day (not including dispersion or free-free 
absorption), for different assumptions on E and tiqbm as labeled. 
The two lower cases are independent of E. 



star creates a low density neighborhood within a denser ISM 
(Whalen et al. 2003). We will refer to the medium relevant 
to shock deceleration as the circumburst medium (CBM), 
and denote its density as Ucbm- 

Considering three different cases, ucbm = 1, 10 -2 and 
(1 + zf, all in units of cm~ 3 , Fig.3 plots the v = 100 MHz 
flux at a fixed time t — 0.01 day (~ 15 min) after the burst, 
not including the effects of dispersion or free-free absorp- 
tion. At this range of v and t, typically the spectra are ris- 
ing with v and the light curves are rising with t. Thus, two 
effects compete in the dependence of the flux on z: time di- 
lation makes one sample the light curve at an earlier time 
with less flux, whereas K-correction makes one observe the 
spectrum at higher frequencies with more flux. This can re- 
sult in a relatively flat behavior if ucbm does not evolve, 
sometimes with a kink when the self-absorption frequency 
is crossed. In constrast, the case of ucbm oc (1 + z) z causes 
such a strong increase in the self-absorption that the flux is 
suppressed rapidly to unobservable levels above z of a few. 
The isotropic equivalent (i.e. beaming-corrected) initial ki- 
netic energy of the blastwave is fiducially E = 10 53 erg from 
the most recent observational determinations (Panaitescu 
& Kumar 2001, 2002, Frail et al. 2001). We also show the 
case of an exceptionally bright GRB with E — 10 54 erg, 
but the difference is apparent only for ucbm = 10 _2 cm~ 3 . 
For either ucbm = lcm~ 3 or (1 + z) 3 cm~ 3 , the flux is al- 
ways in the self-absorbed regime, in which case it is actually 
independent of E. Even in the most optimistic (and rare) 
situation of E — 10 54 erg and ucbm = 10~ 2 cm -3 , the flux 
is at fiJy levels, requiring powerful telescopes for detection. 

Dispersion simply delays the light curve at observer 
time t and frequency v from f(v, t) to f{y, t — At(v)). Con- 
centrating on the most optimistic parameters above, we dis- 
play in Fig.4 the radio light curves at different frequen- 
cies, with and without dispersion for two cases. One is for 
a z — 1 burst with dispersion by an intervening MC of 
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riMC = 10 3 cm~ 3 (which the GRB must reside outside of, 
but near enough so as to be able to sufficiently photoion- 
ize). In this case, free- free absorption by the MC should be 
strong in reality, but we do not account for this here in or- 
der to isolate the consequences of dispersion. The other case 
is for a z = 10 burst with dispersion by the ionized IGM. 
Overlayed is the 3-<r sensitivity of the Square Kilometer Ar- 
ray (SKA) ^, the most sensitive radio telescope currently 
being conceived, assuming that the integration time is 30 % 
of the observed time. 

We see that the effect is to strongly suppress the early, 
power-law rise of the light curves until the delay time At, 
after which the flux rises more sharply than expected with- 
out dispersion. Observationally distinguishing this difference 
is in principle possible, but is likely to be a daunting task. 
The lowest frequency band of the current SKA design is 
at v = 150 MHz, at which At ~ 5 minutes for the MC 
case and ~ 20 minutes for the IGM case. Satellite missions 
such as SWIFT ^ should be able to localize GRB positions 
to arcminute accuracy and distribute them within about a 
minute of the burst. The telescope response time of SKA can 
be short enough to point and begin observing the GRB posi- 
tion within another minute or so. By setting sufficient upper 
limits at the earliest times and characterizing the delayed, 
sharp rise of the light curve, one may be able to determine 
the delay time. The case may be strengthened if this can be 
done at one or more neighboring frequency bands so far as 
At is more than a few minutes. However, establishing the 
absence of a signal near the telescope sensitivity limit is un- 
doubtedly very difficult. The early radio light curve should 
also be significantly modulated by scintillation due to the 
Galactic ISM (Goodman 1997, Waxman, Frail & Kulkarni 
1998), and may be contaminated by emission from the host 
galaxy as well. Most seriously, confusion with other sources 
in the field of view (or even within the angular resolution 
element of SKA; Garrett 2002) is likely to be severe at fj,Jy. 
Hope for distinguishing between different cosmic reioniza- 
tion histories in this way is remote. We conclude that de- 
tecting the effect of dispersion in the radio afterglow emis- 
sion should be extremely challenging but may be marginally 
possible, only for rare, bright GRBs occurring in low density 
environments. 



5.2 Coherent Radio Emission 

The incoherent synchrotron emission of radio afterglows is 
strongly self-absorbed and faint at low frequencies, limiting 
the detectability of the dispersion signature. However the 
observational prospects can be much better if GRBs emit 
strong, coherent radio emission, as predicted by some recent 
models. Very different types of emission from a variety of 
physical processes have been suggested (e.g. Benz & Paesold 
1998, Usov & Katz 2000, Hansen & Lyutikov 2001, Sagiv & 
Waxman 2002). The current observational limits on such 
emission associated with the prompt GRB is quite weak, in 
the 10-100 kjy range, due to the limited sky coverage of 
current facilities and difficulties in filtering out terrestrial 




-2 -1 
log t [day] 



:E=10 54 , n CBM =10" 2 
:z=10, IGM 




-4 -3 -2 -10 12 

log t [day] 

Figure 4. Light curves of radio afterglows for E = 10 54 erg and 
n C BM = 10- 2 cm- 3 , at frequencies v = 10, 30, 100, 300 MHz, 
1, 3 and 10 GHz, from bottom to top, with (solid curves) and 
without (dashed curves) dispersion, a) For the case of intervening 
molecular cloud of tlmc = 10 3 cm -3 at z = 1, neglecting free-free 
absorption, b) For the case of ionized IGM at z = 10. Overlayed 
is the 3 a sensitivity of SKA assuming an integration time 30% of 
t. The dashed portion indicates the range where dispersion effects 
are actually unobservable due to the SKA low frequency limit of 
150 MHz. 



interference (e.g. Dessene et al. 1996, Katz et al. 2003 and 
references therein). 

Particularly interesting for our purposes is the work of 
Sagiv & Waxman (2002). ** Within the framework of the 
standard external shock picture of afterglows, they demon- 
strate that under certain conditions, strong, maser-type co- 
herent emission can emerge on top of the incoherent, self- 



It http://www.skatelescope.org 
I http://swift.gsfc.nasa.gov 



Other models predict emission which is either too faint (Benz 
& Paesold 1998, Hansen & Lyutikov 2001) or at too low frequency 



(Usov & Katz 2000). 
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absorbed synchrotron emission at the onset of shock deceler- 
ation, from either the forward or reverse shock. For plausible 
energy distributions of the injected electrons, this effect can 
be substantial provided that the ratio of the postshock mag- 
netic field energy to electron energy is 1CF 5 , and the CBM 
is sufficiently dense, ticbm ~ 10 4 cm~ 3 . The maser emis- 
sion is characteristically narrow-band around the general- 
ized Razin-Tsytovich frequency, which in the case of the for- 
ward shock in a uniform CBM of ticbm = 10 4 ncsM,4cm~ 3 
is 

v R * ~ 220 li /4 c-L / 1 4 CB^^ 8 n^ Mi4 T^ 8 il MHz , (6) 

where T = 10Tgrb,is is the burst duration, I — 4Lt is a fac- 
tor which depends on the electron distribution, and other 
symbols are as defined in the Appendix (see Sagiv & Wax- 
man 2002 for details) . Although detailed predictions for the 
spectrum and light curve of this component are not yet avail- 
able, the emission should be generally correlated with the 
prompt GRB with durations of ~ 1 minute, and can be as 
bright as ~ 1 Jy for low redshift bursts (E. Waxman, private 
communication) . 

The characteristic frequency of a few 100 MHz in the 
rest-frame is just where the effects of dispersion delay arc 
observationally most interesting for probing the reioniza- 
tion epoch around z ~ 10, which will be observed at few 
10s of MHz (§2; Fig.2). The bright peak flux and fast de- 
caying light curve should allow the dispersion signal to be 
readily distinguished, in the usual manner employed for pul- 
sars (e.g. Cordes & McLaughlin 2003) . Note that this emis- 
sion is stronger for a denser CBM, in constrast to the op- 
posite trend for the self-absorbed incoherent emission. The 
CBM may naturally be as dense as 10 4 cm -3 for the ISM of 
high redshift galaxies or inside molecular cloud cores (§3); 
a stellar-wind CBM is not necessary. Free-free absorption 
by GRB photoionized material may not be a worry here: 
since both the maser emission and UV flash peak at the 
same time, and the forward shock is located in front of the 
reverse shock, the propagation front of the maser emission 
should lead the photoionization front (although the details 
may depend on the differences in the actual light curves). 
The requirement of the low magnetic field may seem more 
stringent but is not precluded, since the mechanism for am- 
plifying the magnetic field above the compressed ISM value 
behind the forward shock is still uncertain. (Values obtained 
from afterglow modeling are in the range es ~ 10 -5 — 10 _1 ; 
Panaitescu & Kumar 2002). In fact, very low magnetic fields 
may be inevitable for the ISM of high redshift (especially 
Population III) galaxies (Abel et al. 2002, Bromm et al. 
2002). All these factors point to this maser emission being 
a valuable probe of the reionization epoch. 

One should first observationally confirm the existence 
of such an emission component from low- z GRBs before it 
can be utilized as a cosmological probe with high-z bursts. 
The best-suited instrument for both objectives may be the 
Low Frequency Array (LOFAR), with its high sensitivity at 
low frequencies, wide field of view, rapid response and effec- 
tive rejection capability of terrestrial signals (Rottgering et 
al. 2002). According to the current design, the sensitivity of 
the whole array at 30 MHz is 68 mjy for an integration time 



of Is and a bandwidth of 4 MHz ft . Assuming that the maser 
emission has a redshifted duration of 10 min (corresponding 
to a rest-frame duration of ~ 1 min at z ~ 10) in the fre- 
quency range IS.v ~ v, and that matched filters accounting 
for the dispersion delay (Cordes & McLaughlin 2003) have a 
detection efficiency of 30%, LOFAR should be able to detect 
this emission at 30 MHz down to ~ 3 mjy, provided that ter- 
restrial and ionospheric interference signals at these fluxes 
can be discriminated. If the mean flux at z ~ 1 is 1 Jy, this 
should be detectable out to z ~ 10 (delayed by 10 hours!). 
Higher sensitivities and limiting redshifts may be achievable 
at higher frequencies as long as that frequency falls in the 
redshifted emission band, and the dispersion delay should be 
measurable as long as the delay time is longer than the red- 
shifted duration of emission. By collecting such information 
for a sufficiently large sample of GRBs at different redshifts 
and determining their redshifts by other methods (Ciardi & 
Loeb 2000, Gou et al. 2003, Meszaros & Rees 2003), one may 
be able to pin down the actual cosmic reionization history 
(Fig.2), which would be a great step forward in elucidating 
the most poorly understood era in the early universe. 



6 CONCLUSIONS AND PROSPECTS 

We summarize the main points of this work. Dispersion 
should be dominated by the ionized IGM for sources at 
redshifts z J> 2, providing a valuable probe of the cosmic 
reionization history at high redshifts, and possibly also the 
low redshift warm-hot IGM. At z 2, dispersion by local 
material ionized by the GRB may be important in molec- 
ular cloud environments, although free-free absorption can 
strongly attenuate the emission at the relevant frequencies. 
For the standard radio afterglow emission, dispersion sup- 
presses the flux in the early-time light curve, the detection 
of which should be very difficult, but perhaps marginally 
possible for rare, bright afterglows with sensitive telescopes 
such as SKA. Dispersion delay should be much more readily 
detectable if GRB emit prompt, coherent radio emission by 
the synchrotron maser mechanism as discussed by Sagiv & 
Waxman (2002). LOFAR may be capable of detecting such 
emission out to z ~ 10 or higher, potentially allowing for 
effective discrimination of different reionization scenarios. 

Dispersion is a unique probe of the IGM in that it is 
sensitive to ionized matter, as opposed to the neutral mat- 
ter probed by most other methods utilizing absorption pro- 
cesses. The main disadvantage is that it can only measure 
the total integrated column density of ionized gas to the 
source, and cannot separate the contributions from differ- 
ent points along the line of sight to the observer. However, 
we have seen that the high-z dispersion signature reflects 
the global IGM and is unlikely be contaminated by other 
local contributions, a crucial asset for probing the cosmic 
reionization history (compare e.g. Barkana & Loeb 2003). 

In this work, we have only addressed the average effects 
of a homogeneous IGM. In reality, the ionized IGM should 
have significant structure, due to 1) density fluctuations re- 
sulting from structure formation, i.e. gravitational evolution, 
which becomes increasingly prevalent at lower redshifts (Cen 

ft http://www.lofar.org 
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& Ostriker 1999, Dave et al. 1999, 2001), and 2) variations in 
the ionization fraction from the inhomogeneous distribution 
of the ionizing sources, which should have a highly compli- 
cated topology before the global reionization redshift (e.g. 
Ciardi et al. 2000). One obvious consequence of the inho- 
mogeneity is to introduce a scatter in the dispersion delay 
around the mean expectation Eq.2 for sources along different 
lines of sight (or a systematic increase for some lines of sight 
with large covering fraction by the warm- hot IGM; §2). Of 
additional interest may be scattering effects: radiation prop- 
agating in an ionized medium with varying density will have 
its path bent due to the varying refractive index (Rybicki 
& Lightman 1979). Fluctuations in the free electron density 
along the propagation path can lead to temporal broaden- 
ing of a transient pulse, as well as angular broadening of 
the image (e.g. Rickett 1990, Cordes & Rickett 1998). The 
magnitude of these effects depends on the spatial distribu- 
tion of the ionized regions and the power spectrum of the 
density fluctuations within them. Estimating their impor- 
tance for propagation through the inhomogeneous IGM is 
nontrivial (e.g. Goodman 1997), and may require compari- 
son with numerical simulation results. If these effects turn 
out to be observationally interesting, further valuable infor- 
mation concerning the density and ionization structure of 
the IGM may be obtainable through the radio emission of 
GRBs. 
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APPENDIX A: RADIO AFTERGLOW 
SPECTRA AND LIGHT CURVES 

We follow standard prescriptions for the synchrotron emis- 
sion from electrons accelerated by a relativistic shock de- 
celerating in the ambient medium (Sari, Piran & Narayan 
1998, Wijers & Galama 1999, Panaitescu & Kumar 2000), 
with particular attention to self-absorption (see also Granot 
& Sari 2002). Aside from the ambient density ncm -3 , here 
taken to be spatially uniform, the principal parameters are 
fiducially taken to be consistent with those obtained from 
broadband modeling (Panaitescu & Kumar 2001, 2002): the 
isotropic equivalent, initial kinetic energy of the blastwave 
is 10 53 i?53 erg, the half angle of jet collimation is 0.1#-i rad 
(Sari, Piran & Halpern 1999; although this is not essential 
for this work), the fractions of postshock energy imparted 
to relativistic electrons and magnetic fields are 0.1e e ,-i and 
0.01eB,-2 respectively, and the spectral index of the electron 
distribution is p = 2.5. We assume an adiabatic shock, for 
which the relation between energy E, bulk Lorentz factor F 
and shock radius r is E = 47rr 2 r 3 nm p c 2 /3, and adopt the 
relation for the observer time t = r(l + z)/4r 2 c. 

The electron distribution is a broken power-law, with 
an injection break at Lorentz factor 7 m = e e Tm p (p — 
2)/m e (p — 1) and a cooling break at Lorentz factor -y c = 
&imi e c/T<T T B 2 t(l + Y), where B = (S2-K€ B m p n) 1/2 Tc. The 
Compton parameter Y is Y = (-^4e e /es + 1 — l)/2 if the 
electrons are fast cooling (-y c < f m ), while it decreases slowly 
with t if the electrons are slow cooling (7™ < j c ) (Sari & 
Esin 2001); we approximate by taking the fast cooling value 
at all times. The power-law index q between these breaks is 
q = 2 for fast cooling and q — p for slow cooling. The ob- 
served characteristic synchrotron frequency associated with 
an electron of Lorentz factor 7 is v — FeBj 2 /ivm e c(l + z), 
and the emitted synchrotron spectrum has breaks at v m 
and v c corresponding to 7 m and 7 C . The flux peaks at 
v v — min(v m ,v c ) and the break frequency above it is 
Vh = max(v m , v c ). The peak flux at v v for a source at lumi- 
nosity distance D is f ViP = (e 3 /3 1/2 m e c 2 )rBr 3 (l + z)/D 2 . 

The evaluation of the synchrotron self-absorption fre- 
quency is based on Granot, Piran & Sari (1999) and 
Panaitescu & Kumar (2000) . Using the self-absorption coef- 
ficient averaged over an isotropic distribution of pitch angles 
(Granot et al. 1999), the self absorption optical depth at v v 



Tp = ipf{p) 



3B7l 



(Al) 



where 7 P = min(7 m , 7 C ), f(p) = (p + 2)(p - l)/(3p + 2) and 
ip = 2 8/ V 5/2 /5r(5/6). This gives, for fast cooling, 

r p = r c = 0.79 x lO^/bXl + Y) 5 e^ 2 _ 2 El A n" 2 , (A2) 
while for slow cooling, 



T P = r m = 3.5 x 10- 7 /(p) ( 2—1 



(A3) 



-5 -1/2 p-l/2„ ( td \ 



where td is the observer time in days. The optical depth at 
an arbitrary frequency is then 



-5/3 



-(«+4)/2 



-(q+4)/2 



V < V p 

v p < v < i4 4 ) 

V h < V 



The frequency at which r v = 1 gives the synchrotron self 
absorption frequency v a , 



3/5 



2/(p+5) / V b 



r b <l< T p 

l-(9+4)/(p+5) (A5) 
T b > 1 



where rt = Tpfab/fp) ^ q+A ^ 2 is the optical depth at Ub- The 
case of Tb > 1 can become important for ambient densities 
n £ 10 4 cm- 3 . 

The spectra f v are broken power-laws with different 
slopes between the different characteristic frequencies. Three 
cases can be distinguished depending on the location of v a 
with respect to v v and for v a < Vp, 
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Va < V < Vp 
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V p < V < Vb 



v b < V 



for v p < v a < v b , 
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V < V p 

V p < V < V a 

(A7) 

V a < V < V b 



V b < V 
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For fast cooling, an additional break may exist at v < v a due 
to the inhomogeneous electron distribution behind the shock 
front (Granot, Piran & Sari 2000), but this is neglected as 
it results in only a small change of the slope. 

The light curves at a fixed frequency v can show several 
breaks of different origin: the passage of some characteristic 
frequency [y a , v v or v^) through v, the transition from fast 
cooling to slow cooling, and the jet transition due to edge 
visibility and sideways expansion. A further break due to 
subrelativistic expansion (Huang, Dai & Lu 1998, Ciardi & 
Loeb 2000) can occur relatively early for a dense circumburst 
medium, but is irrelevant for the present purpose and not 
considered here. 
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